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ABSTRACT 

This work deals with the problem of a boundary layer on a flat plate 
which has a constant velocity opposite in direction to that of the uniform 
mainstream. It has previously been shown that the solution of this boundary 
value problem is crucially dependent on the parameter which is the ratio of 
the velocity of the plate to the velocity of the free stream. In particular, 
it was proved that a solution exists only if this parameter does not exceed a 
certain critical value, and numerical evidence was adduced to show that this 
solution is nonunique. Using Crocco formulation the present work proves this 
nonuniqueness. Also considered are the analyticity of solutions and the 
derivation of upper bounds on the critical value of wall velocity parameter. 
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1 • Introduction 


The boundary layer on the upstream-moving flat plate at zero incidence 
admits of the classical similarity transformation which reduces the relevant 
partial differential equations to the Blasius equation. 

f"' + ff n = 0 
f(0) = 0 

f'(0) = -A, A > 0 
f'(“) « 1, 

where f = \{j(x,y)//(2vx) , \Jj being the dimensional stream function, and v 
the kinematic viscosity, and n = y//(2vx) . This equation can be readily 
integrated once to yield 


i 


• e • , 


f"(n) 


f "(O)exp 


n 

- I 

0 


f (z)dz I , 


f "(h ) = f"(0)exp 


| An i 2 - j f (n - z ) 2 f"(z)dz 


using integration by parts twice. Obviously, the shear stress f"(n) has the 
same sign as the skin-friction at the wall, f"(0). For X = 0, Weyl proved 
the existence and uniqueness using function-theoretical methods. For X 0, 
Callegari and Friedman and Callegari and Nachman found it expedient to work 
with the Crocco formulation, that is, in terms of shear stress g(=f ,r ) as the 
dependent variable and tangential velocity u(=f^) as the independent 


variable: 
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g(u)g"(u) + u = 0, -X < u < 1, 
g'(-A) = 0 
g(l) = 0. 


For X 0, they proved existence, uniqueness, and analyticity of 

solutions to Eq. (2) using an analytical function theory approach. For the 

case X > 0, Hussaini and Lakin proved that a solution exists only for X 

less than a critical value Their numerical results showed nonuniqueness 

for X < X , and the numerical value of X was found to be 0.3541... . In 
— c c 

this work, the nonuniqueness is established rigorously. Also, proof of 
analyticity, and absolute monotonicity etc., is given. Certain analytical 
upper bounds on X are established. 

For convenience, we use the transformation x = u + X to map the 
interval -X < u < 1, to 0 < x < 1 + X. So we consider the equations 


g(x)g" (x) + (x - X) = 0, 0 < x < 1 + X (1.1) 


g'(0) = 0 

g(l + X) = 0. 


( 1 . 2 ) 


2. Analyticity of Solutions 

In this section, the following basic result will be proved: 

THEOREM 1 : There is a range of positive values of X such that the 

positive continuous solution g(x) of the boundary value problem (1.1) and 


(1.2) is analytic on the closed intervel [0,1 + X] 
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This theorem will be proved by considering a sequence of lemmas. The first 
lemma required is: 

LEMMA. 1: The derivative g'(x) vanishes at one and only one point on 
the interval 0 < x < 1 + X. Further , g(x) has its maximum value at this 
point . 


Proof of Lemma 1 : Equation (1.1) can be integrated using the initial 

condition g'(0) =0 to give 


r“ 5 - (2 -» 

Thus, as the initial value a = g(0) > 0, both g(x) and g'(x) are positive 
for 0 < x < X. Also, 

g"(x) = (X - x)/g(x) (2.2) 

is positive for 0 £ x < X and g"(X) = 0. The continuous solution g(x) 

remains positive for X < x < 1 + X, and hence g"(x) is now negative. This 

gives that g'(x) is a monotone decreasing function for x > X. As 
g" ( 1 + X) = -°>, there must thus be at least one point on the interval 
( X , 1 + X) at which g'(x) vanishes. In fact, assuming that g'(x) vanishes 
at more than one point leads to a contradiction, for suppose that g' vanishes 

at both x^ and X£ with x^ < X£. Then, g" would have to vanish at least 

once between these two points which is impossible as g" < 0 for x > X. The 

proof of Lemma 1 is concluded by noting that g"(xj) < 0 implies that 

g(xj) must be the maximum value of g(x). 
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I*EMMA 2: The solution g(x) has a convergent power series expansion on 

the closed interval [x^,l + X]. 


Proof of Lemma 2 : As g(x) is positive and differentiable for 

x ^ 1 + X, equation (2.2) shows that g(x) has derivatives of all 
orders on this interval. Further, expressions for these derivatives may be 
obtained directly from the differential equation (1.1). Induction shows that 
for n _> 1 , derivatives of g(x) satisfy the recursion relation 



(n+1)r g <-» + 1 T [t n ^ 3 ) ♦ («*) 


,001 


g 


( n-k+3 ) 


(2.3) 


where g^) is the k-th derivative of g with respect to x 
the usual combinatorial symbol. 

Let g(x^) = 6, and consider the auxilliary function G(x) 


and (P) is 

V 

defined by 


G(x) = 6 - g(x). 


(2.4) 


Then, as 8 is the maximum value of g(x), G(x) is non-negative for 
Xj < x < 1 + X. Also, for all n > 1, G n (x) = -g( n )(x). Consequently, 
equation (2.1) shows that G'(x) is positive on the interval x _< x < 1 + X. 
From (1.1), 

G (x) = —j— s and G (x) = -j—r 

g(x) g(x) 

are also both positive on this interval. The recursion relation (2.3) thus 
shows that all derivatives of G(x) are non-negative on the closed interval 
[xj,l + X - e] where l+X-x 1 >e>0. Hence, G(x) 


is absolutely 
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monotonic on this closed interval. A theorem of Bernstein [4] now gives 
that G(x) has a convergent Taylor series expansion about the point x 1 
whose radius of convergence is not less than 1 + X - x^ . From the definition 
of G(x), it immediately follows that for |x - xj < 1 + X - x^ g(x) has 
the convergent expansion 


g(x) = l 
n=0 


g (n) (x, ) 

n! 


(X - Xl )“. 


(2.5) 


Application of a Tauberian theorem [5] further shows that the power series 
(2.5) converges at the singular point x = 1 + X to the value g(l + X) = 0 
completing the proof of Lemma 2. 

To establish Theorem 1, it must be shown that for a nontrivial range of 
positive values of X, the power series (2.5) for the solution g(x) of the 
boundary value problem (1.1) and (1.2) converges at the left boundary point 
x = 0. This will be accomplished in Lemma 3. A consequence of this 

convergence will be an expansion for the initial value of g(x) as the series 

» (-I)" *? (n) 

a = 6 + I ^ — - g ( 'Up. (2.6) 

n=2 

tjfmma 3; There exists a positive value X such that if 0 < X < X 
then x^ < (1 + X)/2. 


Lemma 3 gives that the left-hand boundary point x = 0 lies inside the radius 
of convergence of the power series expansion (2.5). Consequently, the 
corresponding solution of the boundary value problem will be analytic. It 
should be noted that the upper bound on x^ given in Lemma 3 is a sufficient, 
but not a necessary, condition for convergence. 



- 6 - 


Proof of Lemma 3: Equation (1.1) may be integrated from 0 to x using 

the identity gg" = (gg')' - (g ') 2 and the initial condition g'(0) =0. A 
second integration from 0 to x^ now gives the result 

*l (x l - 3X) q 2 _ g2 X 1 , 

6 2 + / - Og (OdC. (2.7) 

An upper bound on the right-hand side of (2.7) and a lower bound on the 
maximum point Xj are now required to establish the lemma. 

A lower bound on may be obtained by using (2.1) and the fact that 

g'(Xj) =0 to obtain 



As g(x) is monotone increasing on [ 0^] , g(x) < g(X) on [0,X],but 
S^X) g(x) on [X ,x. ] . Equation (2.8) now gives 


x 1 2X . 


(2.9) 


As g(x) has its only maximum at by Lemma 1, an immediate lower 

bound on g(xj) =8 is 8 > a. A sharper lower bound on 6 can be obtained 
from the expression 


8 



(x 1 - £)(X - 5 ) 
g (?) 


d5 


a + 


~1 

/ 

0 


(X - g ) 2 
g (?) 


d£ 


( 2 . 10 ) 


obtained by integrating (2.1) from 0 to xj. As g(x) < 6 , and by (2.9), 
X 1 * 2 . equation ( 2 . 10 ) now gives the quadratic inequality 
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8 2 - a6 - >_ 0 


which implies 


e > 


a + J a 2 + 8X 2 


/3 


2 2 

A lower bound on B - a which follows from (2.12) is thus 


a 2 2 s 2X 

B - a > — • 


Consider next bounds on the initial value a. Let X = 1 + X. 
integrating (2.1) from 0 to X and using g(X) = 0 gives 


(x - OU - X) 

°*i — tun d5 - 


This relation may be rewritten in terms of strictly positive integrals 


a = / 
X 


(x - 0(5 “ 1) _ r 

g<nD dc A 


(X - 0(x - 5) 

rm 


d? 


which shows 


. r (x - 0(5 ~ O 

“O — gun — dc - 


The convexity of g(x) on [X,X] implies that on this interval 
g(x) >_ g (X)-(X - x). Equation (2.16) now gives that a < (2g(X)) _1 , 
a < g(X), this further implies 

a 2 < V 2 • 


( 2 . 12 ) 

(2.13) 
Then, 

(2.14) 
is 

(2.15) 

(2.16) 
As 


(2.17) 
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Equation (2.15) does not lend itself to the derivation of a lower bound 
2 

on a . However, in the present consideration of analyticity, the required 
bound can be obtained from a relation between a and g which follows from 
the existence proof of Hussaini and Lakin [3]. This proof shows that if X 
is positive and does not exceed a critical value, there is at least one 
initial value a such that a positive continuous solution of the initial 
value problem consisting of (1.1) and the conditions g(0) = a and g'(0) = 0 
exists and has g(X) — 0, i.e., it is a solution of the boundary value 
problem. Further, the solution of the initial value problem will be unique if 

& < 2a. (2.18) 

It must be noted that a unique solution of the initial value problem do«e not 

imply a unique solution of the boundary value problem. This will be shown in 
section 4. 

A lower bound on a 2 follows by using (2.18) in (2.12). The result is 

2 X 3 

0 >.— • (2.19) 

final bound needed for use in equation (2.7) is an upper bound for 
g (x) on the interval [0,xj]. From (2.2), g"(x) is a monotone decreasing 
function on this interval. Further, g M (X) = 0 while the third derivative 
of g is negative when x = X. Thus, g'(x) has its maximum value at x = X. 
This implies that on [0,x^] 

0 < g'(x) < g'(x) = / A—I dc. 

0 su; 


( 2 . 20 ) 



As g(x) > a on [ 0 , A ] , equation (2.20) gives 


0 < g'(x) < ^ . (2.21) 

An upper bound on the integral in equation (2.7) is thus 

X 1 

/ (x - 5 )g' 2 (£)dS < | Xx 2 (2.22) 

0 1 - 8 1 

Use- of (2.13) and (2.22) in equation (2.7) implies 

X l( X l " x ) + 2x3 1 0. (2.23) 

This relation gives that x^ will be less than X/2 for X in the range 
0 < X < X = 0.1176. The sufficient condition for analyticity is thus 
satisfied for a range of positive values of X establishing Lemma 3 and 
Theorem 1 . 

Equation (2.9) implies that x^ cannot be less than X/2 if X > 1/3. 
Indeed, direct numerical solution of the boundary value problem shows that 
< X/2 when X < X = 0.32 and a lies on the upper branch in Figure 1. The 

A 

gap between the values of X and X is associated with fundemental problems 
in obtaining sharper bounds on the initial value a. For example, equation 
(2.15) implies 


x 


1 


< / 
x 


(x - Q(g - X) 
g(?) 


X 

d 5 + / 


U - x) 

g(5) 


X 

d5 - I 

X, 


(5 - X) 2 
g(5) 




a 


(2.24) 
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Individually , the last two integrals in (2.24) are formally infinite, yet they 
must cancel so as to give an order one upper bound. Direct numerical 
calculations show that the upper bound on a 2 is a 2 < 0.219961. The upper 
bound in (2.17) is thus conservative by over a factor of two. 

It must again be noted that < x/2 is only a sufficient condition for 
analyticity. For values of a on the upper branch of Figure 1, solutions of 
the boundary value problem can thus be expected to remain analytic for X 

a 

greater than X. Further insight can be gained by examining parameter values 
for which the condition (2.18), which is sufficient for a unique solution of 
the associated initial value problem, is maintained. Numerical results show 
that (2.18) holds for all values of a on the upper branch of Figure 1. It 
also holds for a on the lower branch of Figure 1 in the relatively small 
range 0.351 < X < X £ and is violated over the remainder of the lower 
branch. The behavior of 8 as a function of a is given in Figure 2. For 
values of X associated with initial values on much of the lower branch of 
Figure 1 , there must thus be serious doubts as to whether solutions of the 
boundary value problem (1.1) and (1.2) are analytic. 


3. An Upper Bound on X £ 

The existence proof of Hussaini and Lakin [3] established the existence 

of solutions of (1.1) and (1.2) for positive values of X less than a 

critical value X . It was shown from (1.1) and (1.2) that X < 1/2. The 
c c 

value of was also determined numerically in that work to be 

X c = 0.3541079... 


(3.1) 



- 11 - 


In this section, additional upper bounds for X £ will be obtained directly 
from (1.1) and (1.2). 

Using the identity that precedes equation (2.7), equation (2.1) can be 
integrated from 0 to x and the result integrated again from 0 to X. 
As g(X) = 0, this gives 

X 2 (X - 3X)_ g a + J ( x - ?)g' 2 (C)dC. (3.2) 

6 1 0 

The right-hand side of (3.2) is intrinsically positive, and thus 


X - 3X > 0. (3.3) 

This relation immediately implies 

X £ V 2 . (3*4) 

To obtain sharper bounds now requires the use of positive lower bounds 
for a 2 and the integral in (3.2). While no additional assumptions are 
required to obtain (3. A), in what follows it will be necessary to assume that 
g < 2a. However, as noted previously, this condition is satisfied on the 
entire upper branch in Figure 1. In particular, it is satisfied in the 
limiting case when X = X £ . 

Let the integral I(x) be defined by 

I(x) = / (X - 5)g' 2 (5)d£ . (3.5) 

0 


Then, as I(X) > 0, equation (3.2) implies 
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X 2 (X - 3X) > 3a 2 . 


(3.6) 


Replacing X by 1 + X and using (2.19) now gives the inequality 


3X 3 + 3X 2 - 1 < 0 


(3.7) 


which yields the improved bound 


X < 0.47533. 


(3.8) 


A slightly sharper bound can be obtained by noting that I(X) > I(X). 
Let 6 = g(X). Then, g(x) <^6 on [0,X], so on this interval 


g' 2 (x) > (2X - x) 2 . 

45 


(3.9) 


This leads to the relation 


I(X) > (5X + 16). 

1206 z 


(3.10) 


An upper bound on 6 now follows from the fact that g(x) > a on [0,X] and 




(3.11) 


In particular, 


j.2 , 2\ J + 1 

5 i 2 


(3.12) 


Use of (3.12) in (3.10) then shows 
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KA) > . >5 »V . 16) ■ 

60(2X + 1) 


(3.13) 


Equation (3.2) now gives 


X 2 (X - 3X) > 3a 2 + 6I( X) 


(3.14) 


which leads to the inequality 


65X 6 + 76X 5 + 10X 3 + 30X 2 - 10 < 0. 


(3.15) 


The solution of (3.15) is 


X < 0.46824 


(3.16) 


which is only a marginal improvement over (3.8). 

Even if the lower bound on I(X) is further sharpened by considering 
this integral on the full interval [0,X], a significant decrease in the bound 
on X is not obtained. Again, this is due to the difficulties associated 
with obtaining sufficiently sharp bounds on the initial value a. 


4. Non-uniqueness of Solutions for 0 < X < X £ 

Using direct numerical results, Hussaini and Lakin [3] have shown that if 
X is positive and less than X^ then solutions of the boundary value 
problem are not unique. For a fixed value of X in this range, as shown in 
Figure 1 there are two initial values a which lead to solutions of the 
boundary value problem. The purpose of this section is to prove this non- 
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uniqueness directly from (1.1) and (1.2). To this end, it is convenient to 
consider the normalized initial value problem 

hh" + t - L = 0, (4.1) 

h(0) = 1, h"(0) = 0 (4.2) 

obtained from the initial value problem for g(x) by taking 

2/1 

g(x) = ah(t) with x = a t. (4.3) 

The parameter L in (4.1) is related to a and X through the expression 

-2/3 

L = a X. (4.4) 


If h(T) = 0 and a(X) is given by 

a = {(1 + X)/T} 3/2 , (4.5) 

then g(X) = 0, so the solution of the initial value problem with initial 
value (4.5) will also be a desired solution of the boundary value problem. 
Equations (4.3) through (4.5) also imply that in terms of T and L 

x - (4,6) 

and 

a = (T - L)~ 3/2 . (4.7) 
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LEMMA 4: Let h^(t) and h 2 (t) be solutions of the initial value 

problem (4,1) and (4.2) corresponding to L values and L£ , 

respectively. Then, if L 2 > L^ , h 2 (t) > h^(t). 

Proof of Lemma 4 : For t « L, h(t) must be of the form 1 + Lt^/2. 

Thus, the lemma holds for small values of t. That it holds for 0 < t T 
can now be shown by contradiction. Let t be the first value of t at which 
h^(t) = 1*2 ("t). As h^ was previously less than t^, this requires 
hj(t) < But, 

L. — t L„ — t L. — t 

h"(t) = - — = - — > ■■■ = h"(t) . (4.8) 

2 h 2 (t) hj(t) hjCt) 

This contradiction establishes Lemma 4. Lemma 4 also shows that if 
h^Ti) = 0 and 1 * 2 ^ 2 ) = 0, then > 0. This implies that: 

COROLLARY: T 2 > Tj . 

The derivative h"(t) is given by an expression analogous to equation 
(2.1). As h(0) is positive,, both h(t) and h'(t) will be positive for 
0 < t < L. This shows that T > L. Consequently, the denominators in (4.6) 
are strictly positive. The following lemma gives a sharper result: 

LEMMA 5: T > 3L. 

Proof of Lemma 5 : Equation (4.1) may be integrated twice from 0 to 


t using (4.2) to give 
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3 2 t 

V 2 h 2 (t) + i- - ~~ = V 2 + / (t - Oh' 2 (Od5. (4.9) 

0 

This implies 

h 2 (t) + I t 2 (t - 3L) > 0. (4.10) 

Setting t = T and h(T) = 0 now establishes the lemma. 

Consider next the behavior of T as a function of L. It has already 
been shown in Lemma 4 that T is a monotone increasing function of L. 

LEMMA 6: T(L) is superlinear in L. 


Proof of Lemma 6 : Let tj be the point at which h'(t x ) = 0. As is the 
case for the original initial value problem in the variable x, there is one 
and only one such point, it lies in the interval L < t < T, and h(tj) is a 
maximum value. 

Equation (4.1) may be multiplied by h' and divided by h to give 


hh" 



0 . 


(4.11) 


Integration from 0 to t produces the result 


V 2 h' 2 + (t - L)lnh(t) - / lnh(5 )d£ = 0. 

0 


Evaluating (4.12) at t^ now shows 


(4.12) 
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One additional bound is needed before demonstrating the superlinear 
behavior of T(L) . The change of concavity of h(t) on the interval [0,^] 
due to the fact that h"(L) = 0 precludes obtaining as a lower bound for h 
on this interval the straight line which passes through the origin and the 
point (tj.Mtj)), i.e., it cannot be shown that h(t) > h.Ct^-t/tj. However, 
for a given L, it is clear that h(t) can be bounded below on this interval 
by a curve of the form 
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h(t )t k 

H(t;k) = (4.18) 

e f 

for a value of k > 1. As k increases, these curves become progressively 
more convex. It should be noted that if H(t,k) provides a lower bound on 
[ °» t : ] for the solution of (4.1) and (4.2) associated with L = U, then, by 
Lemma 4, H(t,k) also provides a lower bound for solutions associated with 
larger values of L. 

This lower bound for h(t) on [0,tj] may be used to obtain an lower 
bound for the integral in equation (4.13). In particular, 

t 

/ lnh(£[)d5 > t. lnh(t. ) - kt . (4.19) 

0 ill 

Equation (4.13) now implies that 

C 1 > £ lnh(t j ) . (4.20) 

Use of (4*17) then gives 

T > t l > 2k ln (4.21) 

The superlinear behavior of T(L) is thus established. 

THEOREM 2: For positive values of X in the range 0 < X < X £ , 

solutions of the boundary value problem (1.1) and (1.2) are not unique. 
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Proof of Therein 2 : Consider the behavior of L as a function of X. By 

(4.4), L(0) = 0. Equation (4.6) and the superlinear behavior of T with 
repect to L shown in Lemma 6 now imply that the graph of X vs L must be 
as in Figure 3. In particular, for a fixed positive X which is less than 
X c , there will be two distinct values of L. By the corollary to Lemma 4, 
each value of L must correspond to a different value of T. Equation (4.5) 
now shows that for the fixed value of X, two distinct values and a ^ 

exist such that the solutions of the initial value problems with these a"s 
are solutions of the boundary value problem (1.1) and (1.2). Solutions of the 
boundary value problem are thus not unique completing the proof of Theorem 2. 
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Figure 1. Values of the parameter a = f ,f (0) for which f (°°) - 


1 as a 


function of X. 
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Values of the maximum value 8 of g(x) as a function of the 
initial value g(0) = a. The dotted line is 8 = 2a. 
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Figure 3. 


The qualitative behavior of the parameter L in the initial 
value problem (4.1) and (4.2) as a function of X. 
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This work deals with the problem of a boundary layer on a flat plate 
which has a constant velocity opposite in direction to that of the uniform 
mainstream. It has previously been shown that the solution of this boundary 
value problem is crucially dependent on the parameter which is the ratio of 
the velocity of the plate to the velocity of the free stream. In particular, 
it was proved that a solution exists only if this parameter does not exceed a 
certain critical value, and numerical evidence was adduced to show that this 
solution is nonunique. Using Crocco formulation the present work proves this 
nonuniqueness. Also considered are the analyticity of solutions and the 
derivation of upper bounds on the critical value of wall velocity parameter. 
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